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ABSTRACT  Presented  are  aspects  of  a  constitutive  model  for  crystalline  metals 
containing  dislocation  and  disclination  defects.  Kinematics  and  balance  laws  are 
developed  simultaneously  at  two  scales  of  observation.  The  macroscopic  kinematic 
description  is  characterized  by  a  multiplicative  decomposition  of  the  deformation 
gradient,  while  the  meso-level  description  follows  from  an  additive  decomposition  of  a 
connection  into  objects  reflecting  dislocations  and  disclinations.  At  the  macroscale,  the 
standard  linear  and  angular  momentum  balances  are  enforced,  while  at  the  mesoscale,  we 
postulate  additional  momentum  equations  for  driving  forces  conjugate  to  defect  densities. 

i 

INTRODUCTION  Experiments  have  indicated  the  formation  of  misoriented  rotational 
and/or  laminar  defect  substructures  within  grains  of  ductile  metals  subjected  to  large 
deformations  [cf.  Ortiz  &  Repetto,  1999].  Following  previous  generalized  continuum 
theories  of  crystal  defects  [Minagawa,  1979;  Pe9herski,  1983;  Bammann,  2001],  we 
develop  a  finite  strain  micropolar  model  wherein  continuously  distributed  dislocations 
and  disclinations  capture  the  physics  of  evolving  defect  patterns  at  multiple  length  scales. 

PROCEDURES,  RESULTS,  AND  DISCUSSION  The  deformation  gradient  FaA  is 
decomposed  multiplicatively  as 

FaA  s  dxa  /  dXA  =  FeaaRa2FpAA  (1) 

where  Cartesian  current  coordinates  xa  -  <f>a(XA,t) ,  and  where  Fe°* ,  Ra2 ,  and  Fp\  are 

the  elastic  deformation,  the  lattice  rotation  due  to  the  time  history  of  disclination  flux,  and 
the  plastic  deformation  attributed  to  the  time  history  of  dislocation  flux.  We  can  further 
decompose  the  elastic  deformation  into  compatible  and  incompatible  parts: 

F‘°, a  =  Fc°aF‘\ ,  2Fc'}\a.h]  =  Fc~'&aj,  -  =  0 ,  (2) 

with  the  subscripted  comma  denoting  partial  coordinate  differentiation.  Disclinations  are 
assumed  to  only  induce  rotation,  i.e.,  RTAa  -  R~ua  [Lardner,  1973;  Pefherski,  1983].  We 
next  introduce  linear  connection  coefficients  and  the  rank  two  metric  tensor  C'«* : 
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(3) 


K  .  F”,F-'-„  +(£,  C*  -  rtV'V 

with  F‘°aF‘-'ac*  the  so-called  “crystal  connection”  of  non-Riemannian  dislocation 
theories,  and  with  ££  a  micropolar  rotation  mode  associated  with  disclinations 
[Minagawa,  1979],  subject  to  the  anti-symmetry  requirements  Qbca=QtcC'da  =  Qb\ay 
The  torsion  and  curvature,  respectively,  of  /”£.  are  written  as  [Schouten,  1954] 

T£  =  F'°a(F>-%F'-t -F‘-'%F'-'\)+ 2Qfa ,  R^  =  -  2 V^+7^ ,  (4) 

where  Va  denotes  covariant  differentiation  with  respect  to  r^.  The  dislocation  density 
a0*  and  the  disclination  density  0‘ * ,  both  referred  to  the  spatial  frame,  are  then  found  as 

r2aah  =  eMT° ,  40*  -  e^e^R^ .  (5) 

Using  Schouten’s  [1954]  identity  and  Bianchi’s  equations  ^.,,=0, 

conservation  laws  for  the  defect  densities  are  written  as  [Minagawa,  1979] 

C‘*ahc'+eaic0*c=  0,  ^=0.  (6) 

We  assume  a  general  free  energy  potential,  per  unit  mass,  with  the  functional  dependency 

^  =  (?) 
where  C'o p  =  FcaagabF‘\  (with  the  Euclidean  spatial  metric  g*)  and  where e =by[ps  is 
a  dimensionless  measure  of  the  statistically-stored  dislocation  line  length  per  unit  current 
volume  ps ,  with  b  the  magnitude  of  the  Burgers  vector.  The  local  energy  balance  and 
entropy  inequality,  under  strictly  isothermal  conditions,  are  written,  respectively,  as 

(lQ^KgL +Z***  +  =  pi*pt,  (g) 

with  p  the  current  mass  density,  e  the  internal  energy,  cr“*  the  symmetric  Cauchy  stress 
satisfying  cr0^  +  pfa  —  pva ,  Xak  and  ^  generalized  meso-level  stresses,  and  -4V  the  Lie 
derivative  with  respect  to  the  velocity  field  Substituting  (7)  into  (8)2  gives 

aai  =  2p{d\p/ dgc,,),  o*g^J\r**r*.-KSZ  0,  (9) 

where  we  define  FaA  =  Ra^FpA.A  and  K  =  pdip/ds,  and  also  provides  the  conjugate 
thermodynamic  force  relations  Xat  ~  pdtf//daai  and  =  pdy/S9  .  We  now  make  the 
more  specific  assumption 

PY  =  PY\\C' op ) + (1  /  2)  ^  (/a2«o4gocgwaa/  +  lt9abgIKgM0cd +cs2),  (10) 
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with  \ffx  the  recoverable  energy  and  fi,  la,  /9,  and  c  311  elastic  shear  modulus,  length 
parameter  associated  with  dislocations,  length  parameter  associated  with  disclinations, 
and  dimensionless  scalar.  Quasi-static  mesoscopic  momentum  balances  are  postulated  as 


[Bammann,  2001] 


****«>  =°» 


g^ac.  +  g^Xlc 


--0. 


(11) 


To  complete  the  kinetic  description,  expressions  for  the  defect  density  fluxes  are  needed: 


1 f'i-f'Hpr}),  ii-ij(M).  MM)-  (12> 

where  {yf},  {/f},  {/l7},  and  jyf]  denote  tensor  functions  of  the  appropriate  rank 

depending  explicitly  upon  the  members  of  the  set  (cr  iXab>£ab’K’X<*.c)  mapped  to 
suitable  configurations.  Spatial  gradients  of  the  dislocation  density  reflecting  the 
contribution  of  pile-ups  to  hardening  are  embodied  by  the  term  %+*  •  Specific  forms  of 
(12)  must  be  developed,  subject  to  constraints  imposed  by  inequality  (%  and  the  force 
balances  (11).  As  an  alternative  to  (12)2,  R\  can  be  found  explicitly  in  terms  of  0 
assuming  stationary  disclinations  [Lardner,  1973],  If  we  further  assume  that  <&  satisfies 
the  integrability  conditions  ^en  the  six  equations  in  (1 1)  are  sufficient  to 

determine  the  six  independent  components  of  Ql  and  (12)3  is  not  needed.  Characteristic 
defect  patterns  are  expected  to  emerge  as  a  result  of  certain  prescriptions  of  la,  lt,  c,  and 
Eqs.  (12)  leading  to  non-convexity  of  the  energy  potential  (10)  [Ortiz  &  Repetto,  1999]. 
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